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We describe Dalitz plots arising in studies of D-meson decays, and what we have learned 
from them in recent years. After some mention of technique and history, we focus on decay 
of the and Df mesons to the if ~if+7r^ final state, and on Df decay to 7r~7r+7r^. Some 
^ mention of other Dalitz plots involving D-meson decays is also made. 
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1 Introduction: General observations on Dalitz plots 



Typically, Dalitz plots are scatterplots of two kinematic variables describing the decay of a pseudo- 
scalar (P) meson to three pseudo-scalars, i.e., P PPP decays. Figures [ij [2) [4) and [5] show 
examples of such plots. Such kinematics should be described by 9 variables (for the three final 
state 3-momenta), but energy-momentum conservation cuts this down to five variables. Further, 
the lack of spin in the parent and daughters leads to a lack of preferred direction in space; thus 
two angles describing the direction of emission of any one daughter in the parent rest frame (RF) 
are inconsequential as is the angle described by the decay plane when rotated with this direction 
as the axis. Thus we reduce the five variables down to two non-trivial variables. The decay 
phase space can be shown to be uniform in either the energies of any two of the daughters in the 
parent RF, or in the mass-squareds of any two daughter pairs. Of course, in each case, the third 
corresponding variable is kinematically related to the two that are chosen for display in the plot. 
Finally, we should note that in case two or more daughters are identical the decay amplitude must 
be symmetric under interchange of those two, which leads to a symmetry in the Dalitz plot which 
can either be explicitly seen or the plot can be folded. 

As we will see below the predicted uniform phase space of Dalitz plots can be exploited to 
bring out the dynamical features of the decays. For instance, resonances in the decay chain stand 
out as peaks in the two-dimensional plot (and are particularly easy to identify and study for the 
mass-squared version). Thus, Dalitz plot studies can be used to find resonances and to measure 
their parameters such as their masses and widths. Finally, since Dalitz plots reveal interference of 
decay amplitudes in such a spectacular way they may be used to search for CP-violating effects 
which also arise due to interference of decay amplitudes. 

2 Fits to Dalitz plots 

The decay amplitude squared leads to the distribution visible in a Dalitz plot. Parameters for such 
an amplitude can be fit using an using either a binned or, preferably, an unbinned fit. 

A typical fit utilizes a model which describes the amplitude as a sum over complex amplitudes 
as follows: 

A^^^Ati (1) 

where ^0 is the S-wave amplitude and each Au is the amplitude describing the i^^ resonance 
of spin I. Such a model is often called the "isobar" model for historical reasons, and the individual 
resonant amplitudes Au are in turn described by expressions such as 

Au = au exp ^ T^u BWu Qu (2) 
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Figure 1: BaBar data [21j: Dalitz plot for Df TT^TT 7T decays. 
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Figure 2: BaBar data [20]: Dalitz plot for D+ K+K' 7r+ decays. 
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where D and TZ^ are form factors of the L)-meson and the resonance i?, and where BW^ and 
Qli are the resonance Breit-Wigner (propagator) and spin factors, respectively. Of course, each 
of these factors is a function of kinematic quantities which are suppressed here in the interest of 
brevity. 

The fits to data require that we integrate the differential decay rate over the phase space, and a 
simplyfying factor is that most of the integrals can be done just once since the complex coefficients 
a^i exp {ij^ii) factor out. Resonance parameters sometimes cause problems in the fitting procedure 
because integrals then have to be done every time these are changed. 

Background contributions to the Dalitz plot have to be estimated either from data or Monte 
Carlo (MC). In principle, the L)-meson mass sidebands in data are a better source of knowledge for 
backgrounds, but do not fit kinematically into exactly the same Dalitz plot. If we apply a i?-mass 
constraint to these background regions to solve the kinematic problem with the plot boundaries, 
the two-body masses are slightly shifted, which causes problems in describing resonances. All these 
effects must be properly accounted for before a successful fit can be done. Variations in detector 
efficiency across the Dalitz plot must also be modeled in some way, typically as a smooth function 
of position in the plot. This is often done using a polynomial or an Artificial Neural Net fit to the 
data. 

Finally, we should elaborate on the 5- wave term ^o- While this is often taken to be just a 
complex constant, more recently experiments such as E791 [7J have experimented with a mass- 
dependent piece- wise linear "function", or set of constants, instead. Surprisingly, this works rea- 
sonably well in the fits, i.e., the resulting constants vary smoothly with mass. This has come to be 
known as the "Model Independent Partial Wave Analysis" (MIPWA) technique. 

3 Early D Dalitz results 

Since there are three flavors of i?-meson (i?^, and D~^)^ and many possible daughter particles 
(neutral and charged tt, K) a large number of Dalitz plots can arise in L)-decays, but only a few 
of these have been studied. Such plots were first made in the 1980's and 1990's by the Mark 
II p.J experiment which observed a "non-uniform" density, by the Mark III [2J experiment which 
published a Dalitz plot analysis, and by the E691 ^ experiment which confirmed the main features 
of the previous fit to K~ti^ti^ decays viz., the presence of a strong non-resonant amplitude, 

and the fact that the isobar model gives a poor overall fit. 

The fit problems were generally traced to the S-wave amplitudes and explored further by the 
E791 experiment. Early work in E791 resolved the fit issues by introducing a /^(800) resonance for 
the K~Ti^ system [4J, and a cr(500) resonance [5j for the 7r~7r+ system. 

Later, following a suggestion by W.M. Dunwoodie of SLAC, E791 pursued a model-independent 
way to compare Dalitz decay data to LASS results on K~ti^ scattering j6j. This E791 study, 
pioneered by B. Meadows, proved to be a success for this MIPWA technique [7j. In the meanwhile. 
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Figure 3: BaBar data: KK S-wave intensity in Dalitz decays as a function of uikk [20]. 



the FOCUS experiment found an acceptable fit to K 7r+7r+ decays using a K-matrix 

description of the S-wave with no cr(500) pole. 



4 Recent D Dalitz Plot results 

More recently, higher statistics data became available from the CLEO, CLEO-c, BaBar and Belle 
experiments. Some of these fits are described below. Also, i?-decay data from these experiments 
have been used to study CP-violation in D decays and also in B decays. 

Shown in Fig. [2] is the K^K~7r^ Dalitz plot from BaBar [20]. Using moments of 

Legendre polynomials of helicity angles, BaBar extracted the S, P magnitudes and relative phase 
in the region (0.99 < mx-K+ < 1.15) GeV/c^ and found that the P-wave is essentially pure 
0(1020), and is therefore described as such to extract a binned S-wave amplitude. Comparing 
5- wave intensities from different Dalitz plots (see Fig. [3]) one finds good agreement except for the 
K^K^ amplitude above 1.1 GeV/c^. Figure [l] shows the BaBar Df 7v^7v~7v^ decays Dalitz 
plot [22. 
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Figure 4: BaBar preliminary data: K^K~ti^ Dalitz plot fit projections assuming no CPV^ 

with the regions used for model-independent comparisons indicated by boxes on the scatterplot. 
The data are represented by points with error bars, the fit results by the histogram. The normalized 
residuals shown below, defined as (N^ata — ^Mc)/ V^MC, lie within a zb5a band in each case, 
with almost all of the values within zbScr. 
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Figure 5: LHCb data: CP asymmetry in progressively smaller bins of the Dalitz plot for 
K^K~7T^ decays. 
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Figure 6: BaBar preliminary data: The difference of the Dahtz plot projections of data (points) 
and the fit (blue band) between the and D~ decays. The width of the band represents the 
±la error expected for our data sample size. 

5 Search for CPV in K^K'tt^ decays 

Most recently there has been an analysis from BaBar on the search for CP-violation in DCS 
K^K~Ti^ decays. Figure [i] shows the combined Dalitz plot for these decays. Searches 
for CP violation {CP\) in charm meson decays provide a probe of physics beyond the Standard 
Model (SM). Singly Cabibbo-suppressed (SCS) decays can exhibit direct CP asymmetries due to 
interference between tree-level transitions and AC = 1 penguin-level transitions if there is both a 
strong and weak phase difference between the two amplitudes. In the SM, the resulting asymmetries 
are suppressed by 0{\VchVuh/VcsVus\) ^ 10~^. Measurement of a larger CP asymmetry could be 
caused by the enhancement of penguin amplitudes in D decays due to final-state interactions [HI |9j 
or by New Physics [TOl [11] . 

The LHCb and CDF collaborations recently reported evidence for the difference in CP asym- 
metry of the time-integrated 7t^7t~ and K^K~ decay rates [I2l EH]. In the decay 
mode K^K~7t^ the CLEO-c collaboration has reported [26j that the asymmetry defined 

by ^ _^ 

AcP = / dmldml I I dmldml (3) 

is found to be (—0.4 ± 2.0^^'^^.^'^)% i.e., consistent with zero. The LHCb collaboration has also 
reported on a search in this decay mode; Fig. [5] shows the asymmetry in various-sized regions of 
the Dalitz plot from that experiment [27]. After taking account of the overall rate, they do not find 
any evidence for CP violation in the differential decay rate. Searches for CP violation in other 
SCS decays with identical transitions c udd and c uss are relevant for understanding the 
origin of CPV [IS [IS]. 

The BaBar analysis does not find any evidence for CP violation in the SCS decay 
K^K-JT^. The integrated CP asymmetry is found to be (0.35 zb 0.30 zb 0.15)%. We find also that 
the asymmetries in four regions of the Dalitz plot are consistent with zero, and that the and 
D~ Dalitz plots are consistent with no CP asymmetry with a probability of 72%, according to 
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Table 1: Preliminary CP violating parameters from the Dalitz plot fit. The first errors are sta- 
tistical, the second are systematic uncertainties which are determined from taking in quadrature 
the errors associated with tracking, the production model correction, the event selection, and the 
Dalitz plot model. 



Resonance r (%) Acf) (°) 



i?*(892)0 


0. (FIXED) 


0. (FIXED) 


i?o*(1430)° 


-dAOtlil ± 4.42 


-G.lltlil ± 1.39 


(p{W20) 


O.Sbtoil ± 0.60 


TAStlfo ± 2.35 


NR 


-14.30l{^|?±5.98 


-2.56l^:?^ ± 8.91 


k(800) 


2.00l|:^^ ± 1-85 


2.10l|^^ ± 1.01 


ao(1450)° 


5.07l:tit ± 9.39 


4.00l|:^^ ± 3.83 



Ax Ay 
/o(980) -0.199tg:i?^ ± 0.084 -0.23ltg:igg ± 0.079 
/o(1370) 0.019tg:g^^ ± 0.022 -0.0045^^;^^^ ± 0.016 



the analysis of the normahzed residuals of the and D~ Dalitz plot divided into 100 equally 
populated bins. Finally, we find no evidence for CP asymmetry in decays through various in- 
termediate states with a study of the two-body mass distributions as seen in Fig. |6) and with a 
parameterization of the Dalitz plot for which the measured CP asymmetries in amplitudes are 
listed in Table [H 
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